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; Abstract 

In the light of recent neutrino oscillation and non-oscillation data, we revisit the phenomenolog- 
ical constraints applicable to three observables sensitive to absolute neutrino masses: The effective 
neutrino mass in single beta decay (rnp); the effective Majorana neutrino mass in neutrinoless 
double beta decay (mpp); and the sum of neutrino masses in cosmology (E). In particular, we 
include the constraints coming from the first Main Injector Neutrino Oscillation Search (MINOS) 
data and from the Wilkinson Microwave Anisotropy Probe (WMAP) three-year (3y) data, as well 
as other relevant cosmological data and priors. We find that the largest neutrino squared mass 



difference is determined with a 15% accuracy (at 2a) after adding MINOS to world data. We also 
find upper bounds on the sum of neutrino masses £ ranging from ~ 2 eV (WMAP-3y data only) to 



' ~ 0.2 eV (all cosmological data) at 2a, in agreement with previous studies. In addition, we discuss 

■ the connection of such bounds with those placed on the matter power spectrum normalization 

parameter erg. We show how the partial degeneracy between £ and as in WMAP-3y data is broken 
by adding further cosmological data, and how the overall preference of such data for relatively high 
values of erg pushes the upper bound of E in the sub-eV range. Finally, for various combination 
of data sets, we revisit the (in) compatibility between current E and man constraints (and claims), 
and derive quantitative predictions for future single and double beta decay experiments. 

PACS numbers: 14.60.Pq, 23.40.-s, 95.35.+d, 98.80.-k 
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I. INTRODUCTION 



One of the greatest challenges in current neutrino physics is to establish the absolute 
masses (mi, m2,m 3 ) of the three neutrino mass eigenstates (z/ 1; z/ 2 , z/ 3 ), for which there is 
ample experimental evidence of mixing with the three neutrino flavor eigenstates (y e , u^, u T ) 
through a unitary matrix U(9i 2 , #23, #13), where 9y are the mixing angles [l|]. 

Atmospheric, solar, reactor, and accelerator neutrino oscillation experiments constrain 
two squared mass differences, 5m 2 and Am 2 (with 5m 2 <C Am 2 ), parametrized as 

. 2 2 2\ mj + m 2 , ( 5m 2 5m 2 2 \ 

(ml, m 2 2 , m 2 3 ) = - + ( - — , +— , ±Am 2 J , (1) 

where the cases +Am 2 and —Am 2 distinguish the so-called normal hierarchy (NH) and 
inverted hierarchy (IH), respectively. The same data also measure 9 12 and 6*23, and place 
upper bounds on 13 (see @, Hf for recent reviews). 

Non-oscillation neutrino data from single (3 decay, from neutrinoless double (3 decay 
(Oz/2/3), and from cosmology, add independent constraints on the absolute neutrino masses, 
through their sensitivity to the observables [l[ 



S 



22 2,22 2,2 2 
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Ci 3 Ci2^i + Ci 3 s 12 m 2 e 92 + s 13 m 3 e 
mi + m 2 + m 3 , 



(2) 

(3) 
(4) 



respectively (with = cos 9ij and Sjj = sin#y, while 02,3 are unknown Majorana phases). 

In a previous work Q , a global analysis of world oscillation and non-oscillation data was 
performed, with the purpose of showing the interplay and the (in) compatibility of different 
data sets in constraining the (m^, m^, S) parameter space. In this work, we revisit such 
constraints in the light of two recent relevant developments: the first results from the Main 
Injector Neutrino Oscillation Search (MINOS) |H[ (see Sec. II), and the three-year (3y) 
results from the Wilkinson Microwave Anisotropy Probe (WMAP) p] (see sec. IV), which 
have a direct impact on Am 2 and S, respectively. In Sec. II we show that, after MINOS, the 
Am 2 parameter is globally determined with an accuracy of 15% at 2a. In Sec. V we discuss 
the degeneracy between S and the cosmological parameter as (which normalizes the matter 
power spectrum 0, 0]), as well as the breaking of such degeneracy with the addition of 
further cosmological data, which strenghten the 2a upper bound on E from ~ 2 eV (WMAP 
3y only) to ~ 0.2 eV (all cosmolo gica l data). Such limits on S are in agreement with 



other recent analyses ja, |7|, la |9|, [10J, [111, |12j . Older (3 and Oz/2/3 results, including the 0z/2/3 
signal claim of jl3j], are also briefly reviewed (see Sec. III). In Sec V we discuss the issues of 
the (in) compatibility among the previous data sets and of their possible combinations and 
constraints in the parameter space {mp,mpp, £). In particular, we show that the case with 
'WMAP 3y only" still allows a global combination (at the 2a level) with the Oz/2/3 signal 
claimed in [l3T |. while the addition of other cosmological data excludes this combination 
(at the same or higher significance level). Quantitative implications for future f3 and 0u2f3 
decay searches are worked out for various combinations of current data (Sees. V and VI). 
Conclusions are given in Sec. VII. 
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II. OSCILLATION PARAMETERS AFTER FIRST MINOS RESULTS 



The oscillation parameters (5m 2 , #12) are essentially determined by the solar and Kam- 
LAND reactor neutrino data combination, which we take from jij (not altered by slightly 
updated Gallium data [14], as we have checked). The corresponding 2a ranges are 



5m 2 = 7.92 (1 ± 0.09) x 10" 5 eV 2 , (5) 
sin 2 fl 12 = 0.314(ll°:g) . (6) 

The parameters (Am 2 , sin 2 #23) are dominated by atmospheric and accelerator v data, 
within the CHOOZ [l5| bounds on sin 2 6*13 (see for details). After the review 0, oscillation 
results have been finalized for the Super-Kamiokande (SK) atmospheric v experiment ^6| 
and for the KEK-to-Kamioka (K2K) accelerator v experiment |T^ |. The SK+K2K analysis 
in 0, however, contains most of the final SK and K2K statistics, and is not updated here. 

New accelerator data in the —>■ disappearance channel have recently been released 
by the MINOS experiment p. The MINOS data significantly help to constrain the Am 2 
parameter [Hf, while they are less sensitive to #23 (as compared with atmospheric data (l6|). 
and are at present basically insensitive to sin 2 #13, although some sensitivity will be gained 
through future searches in the — > u e appearance channel. They are also insensitive to the 
parameters (5m 2 , 612), whose main effect (at #13 ~ 0) is to change the oscillation phase by a 
tiny fractional amount 0.5 5m 2 cos 2#i2 /Am 2 ~ 0.6 x 10~ 2 (see, e.g., [2]), which is negligible 
within current uncertainties. Therefore, the usual two-family approximation — also adopted 
in the official MINOS analysis 2] — is currently appropriate to study MINOS data. 

In the absence of a detailed description of the procedure used by the MINOS collabo- 
ration for the data analysis, we provisionally include their constraints on the parameters 
(Am 2 , sin 2 #23) through an empirical parametrization of the x 2 statistical function, 

2 ( x ~ l \ 2 ( y ■ x °' w - 2 - 74 V 

Xminos - {^) + y ^ ) , (7) 

where x = sin 2 (26 l 23) and y = Am 2 /(10 -3 eV 2 ). This parametrization accurately reproduces 
the official MINOS bounds in the plane (Am 2 , sin 2 2^23) (at least in the relevant region at 
relatively large mixing angles, sin 2 # 23 e [0.3, 0.7]). A more proper analysis will be performed 
when further experimental and statistical details will be made public by MINOS. 

By adding the above x 2 function in the global analysis of neutrino oscillation data per- 
formed in [2j], we obtain the following 2a allowed ranges: 

Am 2 = 2.6 (l±g:if) x 10- 3 eV 2 , (8) 
sin 2 ^ 3 = 0.45(1+°;!) , (9) 

which noticeably improve the previous bounds in 2]. In particular, the 2a error on Am 2 is 
reduced from ~ 24% |2j to ~ 15% after the inclusion of MINOS data. 

The error reduction for Am 2 also reduces the spread of the (dominant) CHOOZ limits 
on sin 2 ^13 [lij], which are Am 2 -dependent. Therefore, MINOS indirectly provides a slight 
improvement of the previous bounds on sin 2 #13 0, which are now updated at 2a as: 

sin 2 ^ 13 = (0.8 +2 : 3 )xl0- 2 . (10) 

Equations (JSJ), ©, ©> © and (fTUj) represent our up-to-date evaluation of the neutrino 
oscillation parameters (at 95% C.L.). 
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3p oscillation parameter constraints 




FIG. 1: Constraints placed by neutrino oscillation data on the parameters (Am 2 , sin 2 #23) sin 2 #13), 
which are affected by the inclusion of the first MINOS results. The results are shown in terms of 
standard deviations from the best fit. Solid (dashed) lines refers to all neutrino oscillation data 
with (without) MINOS. 

Figure H shows the comparison between the limits obtained on (Am 2 , sin 2 #23, sin 2 #13) 
before MINOS (dashed, from 0) and after MINOS (solid, this work), in terms of standard 
deviations from the best fit (i.e., in terms of Ax 2 , where A% 2 is the fitting function to 
all oscillation data [2]). The impact of MINOS on the Am 2 parameter is rather significant, 
and the related uncertainties appear now to scale almost linearly and symmetrically. The 
analogous figure (not shown) for the (5m 2 , sin 2 # 12 ) parameters — which are not affected by 
MINOS data — can be found in 0. 

III. INPUT FROM SINGLE AND DOUBLE BETA DECAY SEARCHES 

Single /3-decay experiments probe the effective electron neutrino mass parameter mp, at 
least in first approximation. The most stringent constraints on rap are placed by the Mainz 
and Troitsk experiments ^8[, the latter being affected by still unexplained time- dependent 
systematics. We assume that both the Mainz and Troitsk results can be taken at face value 
and combined as in 0,0], obtaining at 2a: 

< 1.8 eV . (11) 

If only the Mainz results were used, this limit would be slightly relaxed (rap < 2.2 eV). 
In any case, the above bound is relatively weak and has almost no impact in the following 
analyses, except for one scenario (see Case 1 + in Sec. V A). 
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TABLE I: Our summary of the estimated 0^2/? nuclear matrix elements C mm and their ±2<r errors 
(in logarithmic scale) for seven relevant nuclei, as derived from Ref. |2j|. See the text for details. 



Nucleus 



log 10 (C mm /yr x ) ± 2a 



76 Ge 

82 Se 

100 Mo 

116 Cd 
i28 Te 

i30 Te 
136 Xe 



-13.36 ±0.10 
-12.83 ±0.14 
-13.13 ±0.20 
-12.96 ±0.24 
-14.32 ±0.29 
-12.98 ±0.30 
-13.41 ± 0.56 



Double beta decay searches with no final-state neutrinos have not reported positive sig- 
nals so far 19], except in the most sensitive ( 76 Ge) detector to date (Heidel berg - Moscow 
experiment), where part of the collaboration has claimed a signal at > 4<t level jl3|, recently 
promoted to > 6a level by a pulse-shape analysis [20] . This claim is still considered as con- 
troversial 0] and we shall discuss its implications with due care in the following sections. 
In any case, any claim or limit on the 0u2/3 decay half-life (Ty 2 ) m a candidate nucleus 
constrains the effective Majorana mass (mpp) through the relation 

2 

™h = , (12) 

in the assumption that the GV2/3 process proceeds only through light Majorana neutrinos 
(and not through new interactions or particles) 0]. The relevant nuclear physics is included 
in the matrix element C mm , which must be theoretically calculated. 

Since the above relation is non-linear, and since all the quantities involved (except the 
electron mass m e ) are subject to large uncertainties, we prefer to linearize Eq. (fT2|) and to 
deal with more "tractable" uncertainties by using logarithms as in [il]: 

0u 



In the following, it is understood that linear error propagation is applied to the above "logs" 
(e.g., logC mm ) rather than to the exponentiated quantities (e.g., C mm ). 

In our analysis, we take the theoretical input for logC mm (central values and errors) 
from the quasi-random phase approximation (QRPA) calculations in Ref. |2l|], where it has 
been shown that the nuclear uncertainties can be significantly constrained (and reduced) by 
requiring consistency with independent 2u2j3 decay data (when available). More precisely, 
for a given nucleus, we take the logarithms of the minimum and maximum values of C mm 
from the ±lcr ranges reported in the last column of Table 1 in j^jj; then, by taking half the 
sum and half the difference of these extremal values we get, respectively, our default central 
value and la error for logC mm (2a errors are just doubled). Table shows a summary of 
the log C mm values and their 2a errors derived in this way, and used hereafter. 
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With the theoretical input for C mm ( 76 Ge) from Table HJ the 0^2(3 claim of 
formed in the following 2a range for mpp: 

log 10 (m^/eV) = -0.23 ± 0.14 (0z/2/3 claim accepted) , (14) 

i.e., 0.43 < < 0.81 (at 2a, in eV). See also for our previous estimated range (with 
more conservative theoretical uncertainties). As in 0, we consider the possibility that 
Ty 2 = oo is allowed (i.e., that the claimed Qv2f3 signal is incorrect), in which case the 
experimental lower bound on rapp disappears, and only the upper bound at 2a remains: 

log 10 (m^/eV) = -0.23l^ 14 (0z/2/3 claim not accepted) . (15) 

IV. INPUT FROM COSMOLOGICAL DATA 

The neutrino contribution to the overall energy density of the universe can play a relevant 
role in large scale structure formation and leave key signatures in several cosmological data 
sets. More specifically, neutrinos suppress the growth of fluctuations on scales below the 
horizon when they become non relativistic. Massive neutrinos with S = mi + m 2 + m 3 in 
the (sub)eV range would then produce a significant suppression in the clustering on small 
cosmological scales (see for a recent review). 

The method that we adopt to derive bounds on £ is based on the publicly available 
Markov Chain Monte Carlo (MCMC) package cosmomc [22]. We sample the following set 
of cosmological parameters, adopting flat priors on them: the physical baryon, cold dark 
matter, and massive neutrino densities (u;& = Q^h 2 , uo c = Q c h 2 and Q u h 2 , respectively), the 
ratio of the sound horizon to the angular diameter distance at decoupling (6 S ), the scalar 
spectral index, the overall normalization of the spectrum A at wavenumber k = 0.05 Mpc -1 
and, finally, the optical depth to reionization, r. Furthermore, we consider purely adiabatic 
initial conditions and we impose flatness. 

From a technical viewpoint, we include the WMAP 3y data (temperature and polariza- 
tion) ((J |23| with the routine for computing the likelihood supplied by the WMAP team 
[2^ | . We marginalize over the amplitude of the Sunyaev-Zel'dovich signal, but the effect is 
small: including or excluding such correction change our best fit values for single parameters 
by less than 2% and always well inside the la confidence level. We treat beam errors with 
the highest possible accuracy (see |23j |. Appendix A. 2), using full off-diagonal temperature 
covariance matrix, Gaussian plus lognormal likelihood, and fixed fiducial C/s. The MCMC 
convergence diagnostics is done throught the so-called Gelman and Rubin "variance of chain 
mean" /"mean of chain variances" R ratio statistic for each variable. Our final constraints 
over one parameter (E) or two parameters (£, ag) are obtained after marginalization over all 
the other "nuisance" parameters, again using the programs included in the cosmomc package. 
In addition to Cosmic Microwave Background (CMB) data from WMAP, we also include 
other relevant data semples, according to the following numbered cases: 

1) WMAP only: Only temperature, cross polarization and polarization WMAP 3y data 

are considered, plus a top-hat age prior 10 Gyr < t < 20 Gyr. 

2) WMAP+SDSS: We combine the WMAP data with the the real-space power spectrum 

of galaxies from the Sloan Digital Sky Survey (SDSS) [2^]. We restrict the analysis to 
a range of scales over which the fluctuations are assumed to be in the linear regime 
(technically, k < 0.2/i _1 Mpc) and we marginalize over a bias b considered as an 
additional nuisance parameter. 
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FIG. 2: Constraints placed by different cosmological data sets (1,...,7) on the sum of neutrino 
masses E, in terms of standard deviations from the best fit in each case. 



3) 



4) 



WMAP+SDSS+SN Riess +HST+BBN: We combine the data considered in the pre- 
vious case with Hubble Space Telescope (HST) measurement of the Hubble parameter 
H = 100/i km s" 1 Mpc" 1 @, a Big Bang Nucleosynthesis prior Q b h 2 = 0.020 ±0.002, 
and we finally incorporate the constraints obtained from the supernova (SN-Ia) lumi- 
nosity measurements of [22| by using the so-called GOLD data set. 



CMB+LSS+SNAstier 1 Here we include WMAP data and also consider the small-scale 
CMB measurements of the CBI |28], VSA 0, ACBAR H and BOOMERANG-2k2 
j3l| experiments. In addition to the CMB data, we include the large scale structure 
(LSS) constraints on the real-space power spectrum of galaxies from the SLOAN galaxy 
redshift survey (SDSS) j^J and 2dF survey j^H, as well as the Supernovae Legacy 
Survey data from [33(. 

5) CMB+LSS+SNAsticr+BAO: We include data as in the previous case plus the con- 

straints from the Baryonic Acoustic Oscillations (BAO) detected in the Luminous 
Red Galaxies sample of the SDSS |33 |. 

6) CMB+SDSS+SNAstier+Lyman-a: We include measurements of the small scale pri- 

mordial spectrum from Lyman-alpha (Ly-ct) forest clouds pjjL l36| but we exclude BAO 
constraints. The details of the analysis are the same as those in 0, H?} • 

7) CMB+SDSS+SN A stier+BAO+Lyman-a: We add the BAO measurements to the 

previous dataset. Again, see @, for more details. 
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TABLE II: Input cosmological data sets for seven representative cases considered in this work, 
together with their 2a (95% C.L.) constraints on the sum of neutrino masses S. 



Case Cosmological data set E bound (2a) 

1 WMAP < 2.3 eV 

2 WMAP + SDSS < 1.2 eV 

3 WMAP + SDSS + SN Riess + HST + BBN < 0.78 eV 

4 CMB + LSS + SN As tier < 0.75 eV 

5 CMB + LSS + SN As tier + BAO < 0.58 eV 

6 CMB + LSS + SNAstier + Ly-a < 0.21 eV 

7 CMB + LSS + SNAstier + BAO + Ly-a < 0.17 eV 



The above cases provide a sufficiently rich list of cosmological data combinations, with 
increasingly strong constraints on E. In particular, Fig. |2] shows the constraints on E for 
each case of our analysis, in terms of standard deviations from the best fit of E. None of 
these curves shows evidence for neutrino mass at > la, indicating that current cosmological 
data can only set upper bounds on E. The bounds tend to scale linearly for the richest data 
sets (e.g., for the cases 5, 6, and 7). 

Table II summarizes the bounds on E derived from our analysis in numerical form (at 
the 2a level). Such bounds are in agreement with previous results in similar recent analyses 
which include WMAP and other data, whenever a comparison is possible 0, 0, H, H, EU EH 
lH, and we can derive the following conclusions: 



As already shown in |6( and |9|], the WMAP data alone, in the framework of cosmo- 
logical models we are considering, are able to constrain E < 2 eV at 95% C.L.. This 
limit, which relyes on a single high-quality dataset, should be considered as the most 
conservative. 

Inclusions of galaxy clustering and SN-Ia data can, as already pointed out in the litera- 
ture (see e.g. |4J,|l2( and references therein), further constrain the results. The datasets 
used in the cases 2, 3, 4 and 5 contain different galaxy clustering and supernovae data 
and thus test the impact of possible different systematics. Such cases provide, respec- 
tively, constraints of E < 1.2 eV, E < 0.78 eV, E < 0.75 eV and E < 0.58 eV at 
95% C.L;, respectively. These results are in reasonable agreement with the findings 
of, e.g., p but are slightly weaker than those presented in |8j (case 5 in particular), 
probably as a consequence of the different data analysis method. Just to mention 
few differences, in (if the likelihood analysis is based on a database of models while 
we adopt Markov Chains. Moreover 0] includes variations in the equations of state 
parameter w and running of the spectral index, while here we consider only models 
with w = —1 and no running. Non-standard models with a large negative running 
of the spectral index or with w > — I prefer smaller neutrino masses when compared 
with observations. Finally, [8] uses the same bias parameter for the SLOAN and 2dF 
surveys while the two may be different. 
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• Including SDSS Ly-a data in cases 6 and 7 (as in [2]) greatly improves the constraints 
on E up to S < 0.21eV and E < 0.17eV (95% C.L.). The latter is stronger than the 
one reported in jij], since we are using the updated SDSS Ly-a dataset of p|. This con- 
straint has important consequences for our analyses, especially when compared with 
the mpp claim. However, we remark that the bounds on the linear density fluctua- 
tions obtained from the SDSS Ly-a dataset are derived from measurement of the Ly-a 
flux power spectrum Pp(fc) after a long inversion process, which involves numerical 
simulations and marginalization over the several parameters of the Ly-a model. The 
strongest limit in case 7 should therefore be considered as the less conservative. 

As we can see, cosmology seems to provide the best constraints available on absolute 
neutrino masses. However, sub-eV constraints can be placed only when we consider a com- 
bination of multiple datasets. It is therefore important to check the degree of compatibility 
between the datasets. A possible way to make such check and to better understand, at the 
same time, the data preference for small values of S, is to consider a joint analysis of E and 
of the so-called as parameter, which represents the expected linear root mean square (rms) 
amplitude of matter fluctuations in spheres of radius R = 8/i -1 Mpc. 



Let us briefly remind that the linear (rms) mass fluctuations in spheres of radius R are 
usually expressed throught their power spectrum P(k) in Fourier space (see e.g. [Hf): 



where W(x) = 3(sinx — xcosx)/x 3 is the top-hat window function, and the mass enclosed 
in the sphere is M = 4irpoR 3 /3, with po denoting the background mass density of the 
universe. The matter power spectrum P(k) is fully determined once the cosmological model 
and the corresponding parameters are defined. The effect of massive neutrinos is to reduce 
the amplitude of the power spectrum P(k) on free streaming scales and, therefore, to reduce 
the value of <r 8 . The neutrino free-streaming process introduces indeed an additional length 
scale in the power spectrum related to the median neutrino Fermi-Dirac speed t> me d by 
kj s = AirGp/v^^ oc E/i. For k < kf s , the density perturbation in the neutrinos grows in the 
matter-dominated era while decays for k > kf s . It is possible to show (see e.g. |39j) that for 
wavenumbers smaller than kf s the power specrum is damped by a factor 



where Vt m is the total matter density. Therefore, in general, we expect that the larger E the 
smaller as (or vice versa), namely, that these two parameters are partly degenerate, with 
negative correlations. If og is determined in some way, the degeneracy is "broken" and E 
can be better constrained. 

The value of as can be derived either in an indirect way, i.e., by first determining the 
cosmological model and by then considering the possible values of the allowed P(k), or in 
a more direct way, i.e., by measuring the amount of galaxy clustering on ~ 8 Mpc scales. 
Since CMB observations provide an indirect determination of as, it is important to check if 
this determination is compatible with the other, more direct, measurements made throught 
galaxy clustering. An incompatibility between the two datasets, with, for example, a low 
value preferred by CMB data vs a high value preferred by clustering data, would lead to a 
formally very strong (but practically unreliable) constraint on E (see also 




(16) 



P(fc;E)/P(fc;0) 



-0.087(S/Q m /i 2 ) 



(17) 
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FIG. 3: Joint 2a constraints on the (<7g, E) parameters (95% C.L. for N^p = 1) derived from four 
representative cosmological data sets (1, 2, 5, and 7, as listed in Table II). 

Fig. El shows the results of our joint analysis of the (erg, S) parameters, in terms of 2a 
contours (Ax 2 = 4 for one degree of frreedom) for four representative cosmological data 
sets (1, 2, 5, and 7), chosen among the seven sets studied in this work (see Table II). 
In case 1 (WMAP only), the degeneracy and anticorrelation between the two parameters 
is evident. The addition of galaxy clustering and Ly-a tend to break the degeneracy by 
selecting increasingly smaller ranges for o~ 8 . It turns out that such additional data tend to 
prefer erg in the "higher part" of the range allowed by WMAP only, so that the upper bound 
on £ is pushed to lower values. Should future data hypothetically invert such trend (i.e., 
the current preference for "high" values of as), the upper bounds on E would be somewhat 
relaxed. In any case, the different combinations of current cosmological data in Fig. 3 appear 
to be in relatively good agreement with each other, with significant overlap of the different 
2a regions — a reassuring consistency check. 

V. GLOBAL ANALYSES: RESULTS FOR VARIOUS INPUT DATA SETS 

In this section we discuss the global analysis of all neutrino oscillation and non-oscillation 
data, with particular attention to the (in) compatibility and combination (at 2a level) of the 
various cosmological data sets (numbered as 1, 2, ... 7 in Table |HJ) with the 0i/2/3 signal 
claim of j 1 31 ] - First we examine the most conservative case where a combination is possible 
(case 1 + in Sec. V A), then we discuss the worst case where the combination is forbidden at 
^> 2a (case 7 in Sec. V B), and finally we make an overview of the results for all possible 
intermediate cases considered in this work (Sec. V C). 
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FIG. 4: Superposition of 2a constraints (95% C.L. for ]Vdf=1) placed by /3, Oz/2/3, oscillation, and 
cosmological neutrino data in the three 2-dimensional projections of the (m«, rnpp, E) parameter 
space. Cosmological constraints are labelled as in Table ITT1 for the seven input data sets. The 0v2[3 
lower limit on mpp from the claim in |l3| is indicated as a horizontal dashed line. 



The results discussed in more detail in the following sections can be qualitatively under- 
stood through Fig. EJ which — following the previous work [3] — shows the three orthogonal 
projections of the regions separately allowed at 2a level in the (mp, mpp, E) parameter 
space by: (a) neutrino oscillation data (slanted bands for normal and inverted hierarchy); 
(b) the seven cosmological data sets (numbered lines with "cosmo" label); (c) single j3 decay 
searches (lines with (3 label); and (d) 0u2/3 decay limits (lines with 0^2(5 label). In the 
latter case, the lower limit from the claim in [l3| is shown as a dashed line to remind that, 
formally, it may be accepted or not. In Fig. 01 the various bounds are simply superposed, 
while real combinations of data are performed in the following sections. Note that the global 
combinations can alter the bounds from separate data sets: e.g., the global upper bounds on 
E may be slightly different from those placed by cosmological data only, while lower bounds 
on E (not placed at all by current cosmological data) arise when the oscillation data (i.e., 
the evidence for nonzero neutrino mass) is included. 
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FIG. 5: Comparison (at 2a) between the regions preferred by the 0^2/3 signal claim (horizontal 
band) and by all v oscillation data plus /3 and WMAP-3y data (slanted bands) for normal (NH) 
and inverted (IH) hierarchy, in the plane (mpp, £). The combination of all such data (thick slanted 
"wedge" in the upper right corner of the plot) corresponds to "Case 1 + " in the text. 



A. Case with WMAP-3y data and the 0^2/? claim 

The most conservative choice for the cosmological input is to rely only upon WMAP 3y 
data (case 1 in Table ITT|) . In this case, the upper bound on E is relatively weak (~ 2 eV), 
and the "slanted band" allowed by the combination of oscillation and cosmological data in 
Fig. 4 can reach the horizontal band allowed by the claimed 0z/2/5 signal ^3|, and a global 
combination of all data becomes possible. Figure 5 shows in more detail this situation in 
the relevant plane (mpp, E), where the global combination (from a full x 2 analysis) yields a 
"wedge shaped" allowed region around E ~ 2 eV and rapp ~ few x 10 _1 eV. This case will 
be labelled in the following (meaning: cosmological data set 1, plus 0^2/3 claim). 

For the relatively high v masses implied by the above 1 + case, the Mainz+Troitsk bound 
on mp, included in the fit, slightly tightens the global constraints in the upper part of 
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v oscill. + (3 + 1/2/5 claim + WMAP 3y 




-3 -2 -1 

10 10 10 1 



m,(eV) 

FIG. 6: As in Fig. [51 but in the plane (mpp, mp)- 

the wedge-shaped solution in Fig. 5 (this is the only such case). Not surprisingly, if the 
global combination dubbed were correct, a positive discovery of mp should be "around 
the corner," just below the Mainz+Troitsk bound. For this reason, we find it useful to 
represent in Fig. 6 the same case 1 + , but in the plane {w,pp,mp). A signal in the range 
mp ~ few x lCT 1 eV is clearly expected, and could be found in the next-generation Karksruhe 
Tritium Neutrino Experiment (KATRIN) |4ol |. which should take data in the next decade. 

B. Case with all cosmological data but without the Oz/2/3 claim 

The most "aggressive" limit on £ is placed by all cosmological data (case 7 in Table ITI]) . 
assuming that their combination is correct and that there are no hidden systematics. In 
this case, the upper limit on £ is so tight (~ 0.2 eV), that no combination is possible with 
the 0^2/? claim at face value. The incompatibility is evident in Fig. which shows how 
far are the regions allowed by oscillation and cosmological neutrino data for both normal 
and inverted hierarchy (confined in the lower left corner) and the Oz/2/3 claim (horizontal 
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FIG. 7: Comparison (at 2a level) between the regions allowed by the 0^2/? signal claim (horizontal 
band) and by neutrino oscillation data plus all cosmological data (slanted bands) for normal (NH) 
and inverted (IH) hierarchy, in the plane (m^, £). In this case, (5 decay data constraints are not 
relevant, and oscillation plus cosmological data cannot be combined with the 0^2/3 decay claim. 
This situation corresponds to "Case 7" in the text. 



band). However, one should not be tempted to say that "cosmological data rule out the 0zv2/3 
claim," for at least three reasons: 1) such claim can only be (dis)proved by another 0^2/5 
experiment with greater sensitivity; 2) the Oz/2/3 signal might be due to new physics beyond 
light Majorana neutrinos; 3) independently of the correctness of the Oz/2/3 claim, cosmological 
bounds on neutrino masses should be taken with a grain of salt, given the many systematic 
uncertainties that typically affect astrophysical data, and given our persistent ignorance 
about two main ingredients of the standard cosmological model (dark matter and dark 
energy) . 
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TABLE III: Allowed ranges at 95% C.L. (2a) for the observables (mp, mpp, £), as derived from 
the global analysis of neutrino oscillation and non-oscillation data. The cases labeled from 1 to 7 
(with either normal or inverted hierarchy) include the seven cosmological input data sets discussed 
in Sec. IV, with the decay claim |l3( treated as a one-sided (upper) limit on mpp. Only in 

the case 1 + we fully include the 0^2/3 decay claim in the global combination (the hierarchy being 
irrelevant in such case). See the text for details. 



Case 


TT* V. 

Hierarchy 


mp (eV) 


m m (eV) 


£ (eV) 


1+ 


Any 


[0.42, 0.84] 


[0.42, 0.75] 


[1.3, 2.5] 


1 


Normal 


[0.0057, 0.69] 


[0, 0.66] 


[0.056, 2.1] 


2 




[0.0057, 0.39] 


[0, 0.38] 


[0.056, 1.2] 


3 




[0.0057, 0.25] 


[0, 0.24] 


[0.056, 0.75] 


4 




[0.0057, 0.24] 


[0, 0.23] 


[0.056, 0.72] 


5 




[0.0057, 0.18] 


[0, 0.18] 


[0.056, 0.56] 


6 




[0.0057, 0.063] 


[0, 0.061] 


[0.056, 0.20] 


7 




[0.0057, 0.050] 


[0, 0.048] 


[0.056, 0.17] 


1 


Inverted 


[0.047, 0.69] 


[0.015, 0.66] 


[0.095, 2.1] 


2 




[0.047, 0.39] 


[0.015, 0.38] 


[0.095, 1.2] 


3 




[0.047, 0.25] 


[0.015, 0.24] 


[0.095, 0.75] 


4 




[0.047, 0.24] 


[0.015, 0.23] 


[0.095, 0.72] 


5 




[0.047, 0.18] 


[0.015, 0.18] 


[0.095, 0.56] 


6 




[0.047, 0.075] 


[0.015, 0.073] 


[0.095, 0.20] 


7 




[0.047, 0.064] 


[0.015, 0.062] 


[0.095, 0.17] 



C. Overview of all cases 

From a glance at Fig. it appears that the Oz/2/3 claim is formally compatible (at 2a) 
with oscillation data plus cosmological bounds only if the latter include just WMAP-3y data 
(most conservative case 1). In this case, the global combination dubbed as 1 + in Sec. V A 
becomes possible. In all other combination would be formally possible only by 

stretching (some of) the errors at (much) more than 2a — an attempt not pursued in this 
work. Therefore, except for the case 1 + already examined, we consider the possibility that 
the Oz/2/3 claim is wrong (i.e., that it provides an upper bound only), and combine oscillation 
data with all cosmological datasets (1,...,7) in either normal and inverted hierarchy. 

The results of such exercise are summarized in Table IIII| in terms of 2a allowed ranges 
for the three relevant observables (mp, mpp, E). Such ranges may be useful to gauge the sen- 
sitivity required by future experiments, in order to explore the currently allowed parameter 
space. E.g., the /3-decay experiment KATRIN, with an estimated sensitivity down to ~ 0.2 
eV can fully test case 1 + in Table ITLT1 may test a fraction of the parameter space up 
to case 4, but it cannot access cases 5, 6, and 7, which will necessarily require new /3-decay 
detection techniques to probe the 0(0.1-0.01) eV range for mp. Analogously, in order to 
significantly test the various cases in Table llllj. cosmological data might need to reach a 2a 
accuracy of about ~ 0.1 eV on E (in the global fit), which may be feasible in the future [l2T |: 
on the other hand, 0v2f3 experiments might need to push their sensitivity down to O(0.1) 
eV or even to 0(0.01) eV, which is quite challenging 19]. Cosmological data seems thus to 
be in a relatively good position to run the absolute neutrino mass race. 
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TABLE IV: Allowed ranges at 95% C.L. (2a) for the 0i/2/3 decay half-lives 1% of various candidate 
nuclei, for the same cases as in Table Mil The T®" 2 ranges are derived from the relation Ty 2 = 
m e/( m2 fjf3^'rnm), by using the C mm and mpp input from Tableland Table Hm respectively. The 
shorthand notation X.YeZW means Ty 2 = X.Y x 10 zw yr. For the sake of comparison, in the 
last row we list exis ting experimental lower limits on Ty 2 (at 90% C.L.), as taken from the recent 



compilations in 



3, 



Case 


Hierarchy 


76 Ge 


82 Se 


ioo Mo 


116 Cd 


128 Te 


i30 Te 


136 Xc 


1+ 


Any 


[1.0e25, 3.5e25] [2.9e24, l.le25] 


[5.4e24, 2.3e25] 


[3.4e24, 1.7e25] 


[7.2e25, 4.2e26] 


[3.2e24, 1.9e25] 


[5.2e24, 8.7e25] 


1 


Normal 


[l.le25, oc] 


[3.0e24, oo] 


[5.1e24, oo] 


[3.2e24, oc] 


[6.5e25, oo] 


[2.9e24, oo] 


[4.3e24, oo] 


2 




[3.3e25, oo] 


[8.9e24, oc] 


[1.5e25, oc] 


[9.5e24, oc] 


[1.9e26, oc] 


[8.7e24, oo] 


[1.3e25, oo] 


3 




[8.3e25, oc] 


[2.2e25, oc] 


[3.9e25, oc] 


[2.4e25, oc] 


[4.9e26, oc] 


[2.2e25, oo] 


[3.2e25, oc] 


4 




[9.1e25, oc] 


[2.5e25, oc] 


[4.3e25, oc] 


[2.6e25, oc] 


[5.4e26, oc] 


[2.4e25, oo] 


[3.5e25, oc] 


5 




[1.4e26, oc] 


[3.9e25, oc] 


[6.8e25, oc] 


[4.2e25, oc] 


[8.5e26, oo] 


[3.2e25, oo] 


[5.6e25, oc] 


6 




[1.3e27, oc] 


[3.4e26, oc] 


[5.9e26, oc] 


[3.6e26, oc] 


[7.4e27, oc] 


[3.3e26, oc] 


[4.9e26, oc] 


7 




[2.1e27, oc] 


[5.6e26, oc] 


[9.8e26, oc] 


[8.0e26, oc] 


[1.2e28, oc] 


[5.5e26, oc] 


[6.5e26, oc] 


1 


Inverted 


[l.le25, 3.3e28] [3.0e24, l.le28] 


[5.1e24, 2.5e28] 


[3.2e24, 1.8e28] 


[6.5e25, 4.7e29] 


[2.9e24, 2.2e28] 


[4.3e25, l.le29] 


2 




[3.3e25, 3.3e28] 


8.9e24, l.le28] 


[1.5e25, 2.5e28] 


[9.5e24, 1.8e28] 


[1.9e26, 4.7e29] 


[8.7e24, 2.2e28] 


[1.3e24, l.le29] 


3 




[8.3e25, 3.3e28] [2.2e25, l.le28] 


[3.9e25, 2.5e28] 


[2.4e25, 1.8e28] 


[4.9e26, 4.7e29] 


[2.2e25, 2.2e28] [3.2e25, l.le29] 


4 




[9.1e25, 3.3e28] [2.5e25, l.le28] 


[4.3e25, 2.5e28] [2.6e25, 1.8e28] 


[5.4e26, 4.7e29] 


[2.4e25, 2.2e28] 


[3.5e25, l.le29] 


5 




[1.4e26, 3.3e28] [3.9e25, l.le28] 


[6.8e25, 2.5e28] [4.2e25, 1.8e28] 


[8.5e26, 4.7e29] 


[3.2e25, 2.2e28] [5.6e25, l.le29] 


6 




[9.1e26, 3.3e28] [2.5e26, l.le28] 


[4.3e26, 2.5e28] [2.6e26, 1.8e28] 


[5.4e27, 4.7e29] 


[2.4e26, 2.2e28] 


[3.5e26, l.le29] 


7 




[1.2e27, 3.3e28] [3.1e26, l.le28] 


[5.4e26, 2.5e28] [3.3e26, 1.8e28] 


[6.8e27, 4.7e29] 


[3.0e26, 2.2e28] 


[4.5e26, l.le29] 


Limits (90% C.L.) 


> 1.2e25 


> 2.1e23 


> 5.8e23 


> 1.7e23 


> 7.7e24 


> 2.2e24 


> 4.5e23 



VI. IMPLICATIONS FOR FUTURE 0i/2/3 DECAY SEARCHES 



Concerning Qv2{3 experiments, we have shown in Table IIHI the expected ranges for the 
"canonical" parameter in various cases. However, in practice, the experimentalists 
measure or constrain decay half lives (Ty 2 ) rather than effective Majorana masses (mpp). 
Indeed, some confusion arises in the literature when the sensitivity of different Oz/2/3 experi- 
ments is compared in terms of mpp, but the nuclear matrix elements are not homogeneous, 
or they are not calculated within the same theoretical framework. Therefore, we think it is 
useful to present also more "practical" predictions for half lives T®J 2 in different nuclei, by 



using one and the same theoretical nuclear model [21[ with well-defined uncertainties (i.e., 
by using the matrix elements C mm as summarized in Table UJ) . 

More precisely, the ranges estimated in Table ILTTl together with the theoretical input 
values (and errors) for the nuclear matrix elements C mm reported in Table HI are used to 
estimate the expected 2a ranges for the 0z/2/3 half lives Ty 2 in different nuclei through 
Equation (|13|). As previously remarked, the linear form of this equation makes it easier to 
propagate the uncertainties affecting C mm and mpp. As a result, we obtain in Table II VI 
the 2a ranges for the half-lives T®Y 2 of several candidate nuclei, for all cases (1 + , 1, 2, ... ,7) 
previously considered in Sec. V. 

The last row in Table IIVI shows the current experimental limits on the various half lives, 
as taken from the recent compilations 0, 41 1. It can be appreciated that, except for 76 Ge 



and 130 Te, current limits must be improved by about one order of magnitude (at least) in 
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order to probe any of the cases in Table |l] Even for the most promising nuclei ( 76 Ge and 
130 Te), however, an experimental improvement by a factor of a few in the Ty 2 limits can 
only allow to probe the most optimistic case 1 + , and a fraction of the cases 1 and 2. Probing 
the other cases (3, . . . , 7) will require (much) more than an order-of-magnitude improvement 
with respect to current 0i>2(3 limits. It is thus very important to invest a great effort in future 
neutrinoless double beta decay experiments with increasing sensitivity to longer lifetimes, in 
order to match the typical expectations from current neutrino oscillation and cosmological 
data. 



VII. CONCLUSIONS 



Building on previous work |4j and on updated experimental input, we have revisited the 
phenomenological constraints applicable to three observables sensitive to absolute neutrino 
masses: The effective neutrino mass in single beta decay {mp)\ the effective Majorana 
neutrino mass in neutrinoless double beta decay (jnpp); and the sum of neutrino masses 
in cosmology (£). In particular, we have included the constraints coming from the first 
MINOS results p and from the WMAP-3y data 6], as well as other relevant cosmological 
data and priors. We have found that the largest neutrino mass squared difference Am 2 is 
now determined with a 15% accuracy (at 2a). We have also examined the upper bounds on 
the sum of neutrino masses E, as well as their correlations with the matter power spectrum 
normalization parameter <7g. The bounds range from £ < 2 eV (WMAP-3y data only) to 
£ < 0.2 eV (all cosmological data) at 2a, in agreement with previous studies. Finally, for 
various possible combination of data sets, we have revisited the (in) compatibility between 
current £ and mpp constraints, and have derived quantitative predictions for single and 
double beta-decay observables, which can be useful to evaluate the sensitivity required by 
present and future experiments in order to explore the currently allowed parameter space. 
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